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Abstract We present a class of vorticity functions that will allow for isolated, 
circular vorticity regions in the background of still water preceding the arrival 
of a tsunami wave at the shoreline. 

1 Introduction 

In two seminal papers [TJ [2], Constantin and Johnson proposed a model for 
studying what happens beneath the surface of the ocean before the arrival of a 
destructive tsunami wave at the shoreline. As opposed to other enterprises, 
the fluid does not move here irotationally beneath its surface in a global 
manner, but the water is still with the exception of some isolated, bounded 
regions where it moves with vorticity. Given the particular character of 
the vorticity region discussed in [T], the authors envisaged the possibility of 
more permissive shapes for the boundary of such regions in |2] : the circular 
vorticity region. The analysis has been put on a firm ground via a dynamical 
systems approach in the paper [3]. 

Our intention in this note is to investigate the key features of the technical 
proof from the latter work and, by relaxing them, to propose a new candidate 
for the vorticity function in the circular vorticity region. In a loose manner, 
the new function is a slight deviation from the Constantin- Johnson vorticity. 

2 Phase-plane analysis 

Following j3l Eq. (3.1)], let us introduce the nonlinear ODE system 

£ = -S -/<*>. r>0 ' (1) 



1 
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where / : R — > K is a continuous (vorticity) function with some extra re- 
strictions to be described in the following. First of all, / has three zeros: 0, 
u G (0, 1] and u\ = —u$. See Figure 1. 

The energy and its decay. We define the "energy" E(ip,f3) = 4r + 

Jo f{u)du and ask also that / be odd, meaning that E(±tp, ±/3) = E(ip,/3), 
together with lim f{u) = ±oo. The latter restriction leads to 

u— >±oo 

i r>> 

lim - / f(u)du = +oo (2) 
by means of the L'Hopital rule. 



Figure 1 

If E(r) = E(ip(r), (3(r)) for some solution of the system ([1]), with 

r G (0, Too) and r ro < +oo, then 

B 2 

E'(r) = pp' + f(M' = -—, r>0. 

r 

This means that the energy E is monotone non- increasing along the solution 
(■0, B). Obviously, for the three equilibria (0, 0), (uq, 0), (u\, 0), we have con- 
stant energy throughout (0, +oo): for the null solution, and f Q u ° f{u)du = 
Jo" 1 f( u )du for the other two. To ensure that the energy E{r) is strictly de- 
creasing along any non- equilibrium solution whenever it is positive, 
one way is to impose that f Q Ul f{u)du < 0, where i G {0, 1}. This (strict) 
monotonicity restriction is a key feature of the estimate [3] Eq. (3.13)]. In 
our case, following Constantin's technique, if E{ri) = E{r 2 ) > 0, which 
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reads as = E(r 2 ) — E{ri) = — ^rdr for some T\ < r 2 from (0, r^), then 
(3 = in [rx,r 2 ]. Taking into account the second of equations (jTJ, we get 
f{ip{r)) = throughout [ri,^] which means that either if) — Uq or if> = ui 
in [rx,r2]. However, for any of these solutions we have E(r{) = E{r-i) < 0, 
which is a contradiction. 

Another way of establishing this decay of the energy along non-trivial 
solutions is to ask that / be C 1 in R — {0}, as in the case of Constantin's 
vorticity function, namely 




(3) 



Now, given (ip, (3) a non-trivial and non-equilibrium solution, that is ip 2 +j3 2 > 
everywhere in (0^^) and tp{r) £ {0,m ,Mi} when /3(r) = 0, we have 



2/3/3' /3 2 
E = h — = - 



— + / m 

r 



(4) 



and 



E" 



2/3 



(5) 



Obviously, /3" exists for all the values r > with ip(r) ^ 0. Since the solution 
(if), /3) is continuous, then in the (small) vicinity of any r + > such that 
/3(r + ) = 0, meaning ip(r + ) {0,u ,Mi}, the function ip has constant, non- 
null, sign. So, by means of © , © , we get E'{r + ) = E"(r+) = 0, E'"(r + ) < 0. 
The Taylor expansion E{r++h) = E{r + )+E'{r + )h+E"{r + )%+E'"{r + )^ + 
o(\h\ 3 ) yields E(r + + h) < E(r+) when h 3 = and h > 0. If /3(r) 7^ then, 
obviously, the energy E decreases strictly on a small right-neighborhood of 
r. 

Local existence I. Assume that the continuous function / from ([T]) 
satisfies the following restriction: for any a > 1 there exists 77 = 77(a) G (3, ~1 
such that 



1/(01 — V a when £ 6 



1 - T. a, 



1 + f M 



(6) 



In particular, for the vorticity function from ([3]) we have |/(0I < 10 + 
(l + j + y/1 + I) a <na whenever a > 1 and 77 > ^. 
We claim that the integral equation 



10 < 



^(r) 



Tf(lJ>(T))dTdZ, r > 0, 



(7) 
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has a solution in X = C([0, 1], M). The proof of our claim is based on the 
Schauder fixed point theorem and on the Ascoli-Arzela relative compact- 
ness criterion. Notice first that, by means of the L'Hopital rule, we have 
lim | J t } '(if)(r))dr = lim[£/ (■?/;(£))] = 0, which means that the function 

£ !->• I J rf(if>(r)) can be prolongated backwards to zero as a continuous 
function. 

Now, consider D the closed ball of radius ^ and center a of the Banach 
space B = (X, \\ ■ W^). Obviously, the ball is convex. Let T : D — > X be an 
operator with T(if)) given by the right-hand member of the equation (J2J) for 
any if) G D. We have the estimates: (i) boundedness ofT(D), namely 

\T(if))(r) -a\< f\fr- ^adrdi = ^ • r < if; G D. 

Remark that we have obtained actually that T(D) C D. (ii) Equicontinuity 
ofT(D), which follows from 

IM'(r)|<- [ r r-r,adT<^, re (0,1]. 
r Jo 1 

According to the Ascoli-Arzela criterion, the set T(D) is relatively compact 
in B, yielding that the operator T transports bounded sets into relatively 
compact ones. 

Finally, (Hi) the operator T is continuous. Since the function 



./ 



is uniformly continuous, for every e > there exists 5 = 5(e) > such that, 
given Vi, ^2 e D with - ip 2 \\oo < S, one has \f(if>i(r)) - f(if>2(r))\ < e, 
r G [0, 1]. We have 

Wi)(r) - T(if; 2 )(r)\ < jT i jf r • e drd£ < |, r G [0, 1]. 

According to the Schauder fixed point theorem, the operator T has a 
fixed point if) € D. In particular, if>(0) = a. Also, ^ G C 1 ((0, 1],R) and 
-0'(r) = - JT rf(ip(r))dr, r > 0. By means of L'Hopital rule, lim^r) = 0. 

This means that if) G C 1 ([0, 1],R) verifies (j7j) and, given its smoothness, 
yields the solution (if), if/) G C([0, 1],M 2 ) of the system (pQ) with the starting 
point (^(0), ^'(0)) = (a, 0). In particular, by being an element of D, the 
function if) satisfies the inequality 

i>{r)>(l-^)a>l (8) 
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Local existence II. We are interested here in the backward existence 
and uniqueness of the solution to (jTJ). Assume that, in addition to (JBJ, 
there exists L = L(a) G (0, §) such that - /(£ 2 )| < - &| for 

all £1, £ 2 ^ [(l — f) a , (l + 4) a ]- m the particular case of ([3]), we have 

l/'K)l<i + ^<i + 7 ^<i + ^ sl + ^ < i- 

Take T > 6 and introduce the system of integral equations 



tjj(r) =ipr- J* P(s)ds, 



r<T, ip T , f3 T G R. (9) 



Here, i/j t = a, \Pt\ < f • For X = C([VT 2 - 1,T],M 2 ), let D = {(V>,/9) e 
X : |V(r) - a| < f , |/3(r)| < 2|/3 T | + rja, r G [VT 2 - 1,T]}. Obviously, £> 
is a closed, convex, with non-void interior subset of B = (X, \\ ■ \\oo), where 

||(^^)||oo = niax{||V||oo,||/3||oo}. 
Set A G (1, 3) and k > such that 



(T + VT 2 - l)lnA (>ll-lnA) 

hlA > fc > max {AlnA,L ( ^AlnA + \ ) } (10) 



T - y/T 2 - 1 L \4 2 

and introduce the metric 

d((V>l,/3l),(V>2,/3 2 )) 

sup [max{|^i(r) -^ 2 (r)|,|/3!(r) -/3 2 (r)|}-e- fcr ] . 
re[v / r 3 =T,T] 

Notice the by-product inequality 11 • In A > L (| • 3 In A + |). This is equiv- 
alent to if^fz2 > Since the function x i-> 1 ^T 2 is increasing in (0,ln3) 
and i5^fz2 > § > -^( a )> suc h a constant A exists always. 

Then, .M = (D,d) constitutes a complete metric space. Introduce further 
the operator T : D — >■ X with T(?/>, /3) given by the right-hand member of 
the system © for any /?) G D: T(if>,P) = (Tx(ip, /3), T 2 (ip, /?)). 

We have the estimates 

T T r T 



W,P){r)\ < \Pt\- -^== + vw — [ j r Vads 

< 2\fr\ + 2 V a(T - VT 2 ^!) = 2\Pr\ + - 2 f__ 

i + V-t — 1 

< 2|/3 T | +r/a 
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and respectively 



|.7 ! ,(r. >')(/') -u\ < C {2\p T \+r 1 a)ds< 2| '''' '" ''" 



T + y/T 2 - 1 



2 • f + r]a _ na 
< 5 ~ T' 



They read as T(D) C D. 

Given (V>i, A) £ D, i E {1,2}, we get 

ITx^^OM-T!^,^)^)! < / e fcB ds-d((^i,/9i),(^^2)) 

k(T—r) i 

Recall the elementary inequality e x < 1 + Ax for any A > 1 and x G [0, In A]. 
We have x = fc(T — r) < k(T — \JT 2 — 1) < In A by means of the first of 
inequalities (ITDjl. So, e fc(T ~ r ) — 1 < AA;(T — r) and we get 

iTx^i^OW-Tx^^WI-e-^ < A(r-r).d((^i,/3i),(^,/32)) 



where r £ [\/T 2 — 1, T\. 
Further, 



L ' T 



r J r 

and, since 

1 r T kTJ Te kT -re kr 1 f T krJ Te kT - re kr 
re dr = : — / e dr < 



kr kr J r ~~ kr 

< kr _ Te fe ' T - r ' - VW^l kr _ Te k ^ — (T — 1) 

~ e Ayr 2 - i ' ( ' avt 2 - i 

< ckr T-Xk(T-r) + l <ekr T-AlnA+1 

*VT 2 - 1 fcVT 2 - 1 

efcr _|yr 2 ^T.AlnA + ^i_^ |AlnA + i 



fcVT 2 - 1 k 
we conclude that 

tZ(r(Vi, fa),T(1> 2 , fa)) < C ■ ft), (V>2, ft)), 
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where 



_ max {A In A, L (SAJsA + 1)} 
^ ~ k < 1 



by means of the second of inequalities fflO|) . Being a contraction, the operator 
T : M. — > M. has a unique fixed point (tp, (3) G D which is the solution of 
system Q in [\/T 2 — 1, T]. See Figure 2. 



Figure 2 

In the particular case when ip T = a = u = 1 and (3t = 0, see e.g., 
the vorticity function ([3]), the uniqueness of solution leads to ip = M o m 
[\/T 2 — 1, T]. This remark has an essential consequence: if there exists a 
non- equilibrium solution (ip, (3) of system that will eventually reach the 
equilibrium (uq, 0) then either the equilibrium is attained in r < 6 "units of 
time" or it is attained in "infinite time", that is lim (ip(r), (3(r)) = (w ,0). 

r— >+oo 

In fact, let us assume that [^(j")] 2 + [/3(V)] 2 > u$ in [0,6] and there is some 
T > 6 with ip(T ) = u and (3(T ) = 0. Introduce T_ = inf{T > 6 : ip(T) = 
U0) p(T) = 0}. Then, T_ > 6 and ^ (T_) = u , /9(T_) = 0. Thus, the 
solution (■0, (3) is constant in [^T 2 — 1, T_], and so we have ^(a/T 2 - 1) = 
Uo, /3(a/T 2 — 1) = 0, which contradicts the definition of TL. 

Polar coordinates (the Priifer transform) I. Let (ip,/3) be a solu- 
tion of the system (TTJ) starting from (a, 0), with a > 1. The solution exists 
throughout [0, 1] and we have the estimate ^(1) > | by means of (jSJ). Since 
[^(l)] 2 + [/?(1)] 2 > 0, we introduce the new variables 

<4) = Rcos9, (3 = Rsin6, R = R(r), 9 = G(r), 

which make sense on a small right-neighborhood of r = 1 via the Peano 
existence theorem for (JTJ. Since E(r) < E(l) for as long as the solution (?/>, 0) 
exists to the right of 1 (recall that E'{r) < 0), the estimate (j2J) shows that 
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both of the quantities ip(r), /3(r) are bounded on their maximal interval of 
existence. According to Wintner's non-local existence theorem, this implies 
the global existence in the future for the solution (i/j,{3). 
In polar coordinates, the system (0Q) reads as 



1 - f • sin 29 + 



2r 

R = -&- sin 2 9+ [Rcos9 - f(Rcos9)} sin 9, 



cos 2 9 - /(fic °y )cos ' 



r>l. (11) 



We are interested here in generalizing [31 Eq. (3.13)]. To this end, assume 
that there exists an odd continuous function g : R — > R such that f(tp) = 
ifj — g(ip) and Jq g(u)du > , ^6l, for some A G (|, l). In particular, 



for 03]), g(ip) = y\ip\ ■ sgntp and A = |. Also, for the vorticity function 
/(V') = — \i[>\ a • sgnip, where a G (0, 1), we have A 
The first of equations (ITT)) reads as 



l + Q 

2 



„. sin 2$ il>q(il>) , s 

? = r>l. (12) 

Assume that i£(r) = E(if)(r), /3(r)) > for r G [l,r ++ ), where r ++ < 
+oo, and ipg{ip) > for all ip G R. Then, we have the estimates 

< £/(r) = / g(u)du < — 



o 



2A 



which, via (fT2"j) . lead to 



-l-l<0'(r)<-(l-A) + ^ re[l,r ++ ). (13) 

To add some geometric flavor to this analysis, notice that the "problem- 
atic" term in the first of equations (ITTj) . that is , reads for fl3j) as ■ co ^ — , 
see [3], Eq. (3.11)]. Let (i(r)^(r), t(r)/5(r)) be the point of intersection be- 
tween the algebraic curve E(if), 0) — and the straight line passing through 
the origin O of the phase plane Otpf3 and the current point of the flow (so- 
lution) (ip(r), /3(r)), r G [l,r ++ ). See Figure 1. Thus, u — y - and so we 

i /TT Icos0| 3 / 2 

have Vkl = ^7^- 

Polar coordinates II. Let us impose a new restriction of the last term 
of (TT2]), that is 

< imn. < U < r 14 

if _ R 2 ~ R v ' m _ ' v ; 

for some c G [0, 1), f G (0, 1). In the particular case of (JHJ), we have c = 
and v = h- 



Nonlinear ODE's in tsunami modeling 



9 



Take 5 > e > (1 + c) l ' v - 1 > and a > 8(1 + 5). We are interested here 
in the eventual encounter of the trajectory {(ijj(r), /3(r)) : r > 0} starting 
from (a, 0) with the ring l + e<R<l + 5, see Figure 1. According to ([8]), 
this encounter can happen only for r > 1, so we can use the equation (fl2|) to 
investigate it. We claim that liminf R(r) < (1 + c) 1 /", meaning i/ie solution 

will reach the ring in "finite time". 

For the sake of contradiction, suppose that liminf R(r) > 1 + e. Then, 

r— >+oo 

by means of ffl2l . we have 

2 (1 + e) v ~ 2r R v ~ 

1 1 + c 
< -H h t r-i r > r_ > 1. 

2r_ ~ ~ 

Obviously, -R(r) > 1 + e for all r > r_. Since (1 + c)(l + < 1, we can 
take r_ greai enough to have —1 + ^ — I- (1 + £)~ u = —f] < 0. By now, 
~~ I ~~ (i+e) v —0 — ^ throughout [r_, +oo). 

Recalling fl2]), observe that the quantity E(r) is bounded from below 
and monotone non-increasing. This yields lim E(r) G R, and so we have 

r— >+oo 
/3 2 

latter improper integral leads to 



lim E{r) — E{rJ) = f °° E'{r)dr = — f °° —dr. The convergence of the 

r— s>+oo r ~ r— r 



N2 f + °° sin 2 6 , f + °° (5 2 , 
0<(l + e) 2 / rfr < / — rfr < +oo. 

The double inequality — |(r — r_) < 0(r) — #(r_) < — w(r — r_), r > r_, 
shows that lim 9(r) = — oo. The function 9 being strictly decreasing, we 

r— >+oo 

deduce the existence and uniqueness of the increasing, unbounded from above 
sequences (r+) n > no , (r~) n > no , with r+ > r~, given by the formulas 

7T 3"7T 

#( r n) = T - 2mr + 6 '( r -)' ^( r n) = -j-2mr + 0{rJ), n > n , 



where no > 1 + 



\8(r-)\ 
2tt 



Remark also that 



0(rti) = \ ~ 2(n - 1)tt + 0(r_) = 0(r") + y > e(r"), 

and, accordingly, r^_ : < r~, n > no + 1. 
We have the estimates 

— dr > W ^dr > V ^— pL. (15) 

ri=no " n=no 



Nonlinear ODE's in tsunami modeling 



10 



Taking into account that 

c 



3 

2 + (1 + e)" 



(r+-r-) < e'(Q(r+-r-) 
= 6{ri)-6{r-) = 



7T 

2 



< -r/(r+-r n ), £ n e(r n ,r+), 



which yield 



7T 

2?7 



> r — r > 

— n n — 



7T 



2c ' 



n > no, 



we deduce that 



(l+e)- 



(16) 



no 



+ Y ( r n~ r n-l)+ Y ( r n~ r n) 



n=no+l 



> r, 



«■<) 



Y ( r n ~ r n)> r n 



71 



n=ng+l 



3- 



(l+e)" 



for any m > n + 1. Notice the by-product inequalities 2w < 2 < | 
3 — tt- < 3 — Trrru- Further, we have 

1+c ' 



3+c 
+c 



7T 



2nvr = 0(r+)-0(r_) = / (// n )(r+-r_) 
< -7/(r+-r_), u n G(r_,r+), 



which yield 



2n?r - ? 



V 



+ r_ > rj~, n > uq. 



(17) 



In conclusion, via (|T5|) - (|T7|) . we obtain that 

+0O 



^oo •„ 2 



sin^ 



4rj 



n=?io 



7T?7 



E 

n=no 
+oo 



3- 



4 [ ^ + r_ 



E 



12 - 2ira + (T + *■-) 



+ OO. 



We have reached a contradiction. 

Polar coordinates III. Escaping E = 0. The situation depicted in 
Figure 3 below is inspired by the case of (JHJ). There, see [3], in quadrant I, that 
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is when 9 G (0, |)+2tt-Z, the algebraic curve E(ip, (3) = /3 2 -||V'| 3/2 +V' 2 = 
is concave, the point (0, 0) is singular and at the "smooth peek" (ip + , 0), where 
ip + = y, the curvature is |. The energy of the non-null equilibria is also 
negative, J " ! f{u)du < 0, where i G {0, 1}. A consequence of the concavity 
reads as follows: in quadrant I, each oblique straight line splits the region 
E < into two disjoint parts, the "+ " and the "— ", see Figures 1, 3. 



Figure 3 

A solution (ip, (3) of ([I]) starting from (a, 0) for some very great a > 1 
will reach the ring 1 + e < R < 1 + 5 in finite time r. Since ip + > 1 > m , 
the positive quantities e, 5 can be taken small enough, meaning such that 
ip + > 1 + 5 > 1 + e > (1 + c) 1/u > 1, for the ring to intersect the region 
E < of the phase plane as presented in Figure 3. The decay of E(r), 
detailed previously, has the following consequences: if the point (ip(T), f3(T)), 
for some great T > 1, is on the algebraic curve E = then the trajectory 
will cross (transver sally) from outside to inside the algebraic curve; if the 
trajectory will intersect the algebraic curve in its (smooth) peek (ip + , 0) then, 
again, the trajectory will enter the right "lobe" of region E < 0. 

Suppose, for the sake of contradiction, that the solution (ip, ft) remains 
for ever outside the region E < 0, that is liminf E(r) > 0. As in Figure 3, 

i — >+oo 

this means its trajectory will intersect the horizontal axis outside the interval 
The estimates ([TBI . (P3J, namely 



show that the trajectory will rotate around the origin O infinitely many times 
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and it will also intersect the ring infinitely many times. Taking into account 
the possible positions of the points from the trajectory that lie inside the ring, 
we have four (in symmetrical positions) "siblings" of the situation depicted in 
quadrant I. There exist, accordingly, two increasing, unbounded from above 
sequences (r+) n > no , (r~) n >„ , with r+ > r~, given by the formulas 

0(r+) = #1 -2rar + 0(r_), 6(r~) = 9 - 2nn + 0(r_), n > n , 

for some great integer tiq. As before, since 2it — (6q — 9{) > 0, we have 
r^_! < r~ when n > no + 1. Notice also that R(r) > 1 + e > 1 whenever 
re[r-,r+], see © Eq. (3.18)]. 

When the solution escapes from quadrant I while remaining outside the 
region E < 0, its possible trajectories in quadrant IV read as in Figure 4. 
Repeating the argumentation from Polar coordinates II, we reach a contra- 
diction. So, lim E(r) < 0, which means that the solution (ip,/3) will enter 

the region E < in "finite time". 



Figure 4 

Let us return to the estimates (TT3"]) . (|T4|) . Introduce C G (0, 1 — A) such 
that r ++ > 2jYi~\wi > 1. If a solution (ip, (3) is still outside the region E < 

for large values of r, meaning for r > 2 [(i-a)-c] ' ^ nen we nave the estimates 



3 



2 (1+e) 1 - 



< 6' < — ( throughout 



2[(i- 1 a)- C] » r ++ ) • The "+/-" splitting of 
the region E 1 < displayed in Figure 3 implies that, if the solution (ip, (3) will 
not leave quadrant I ever again, then its trajectory will cross from outside 
to inside the algebraic curve E = but it will not encounter the origin O as 
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the angle 9 is decreasing with non-null "angular velocity" 9' ! In conclusion, 
if the solution (if), 0) has not reached the origin O in a short period of the 
time r then the only possible "entry windows" are quadrants II and IV. See 
Figure 4. 

Transversality of flow on the vertical axis. Assume again that the 
energy of the non-null equilibria is negative, J™ 1 f{u)du < 0, where i G {0, 1}, 
and the algebraic curve E = intersects the vertical axis [ip = 0] only in the 
singular point O. If the solution (tp, 0) of ([I]) is outside the region E < 
then we have 

V/ = /3^0, 0? = -K r>l, 

r 

whenever if)(r) = 0. Take > 1 such a zero of the function ip. Applying the 
local inversion theorem, we deduce the existence of a C 1 -function r = r{ip) 
with r : (— £*, £*) — >■ (r* — 5*, r* + <5*) C (1, +oo) for some small £*, 5* > 0. We 
have now a new ODE for characterizing the other function = /3(r(/0)) = 
(3(i/j), namely 

/3(0) = & ^ 0. 

The function /3 i— >■ 4 being locally Lipschitzian in K — {0}, the solution 
of this new ODE exhibits continuous dependence of the data. The same 
issue has been dealt with in [3J Eq. (3.12)] by means of polar coordinates. 
The topological arguments from [31 Section 3] can now be applied verbatim 
to establish that there exists some a > 1 great enough in order to have a 
solution (i/),f3) that will reach the null equilibrium O in finite time r. 



3 An example of vorticity function 

Introduce the function 



f(u) = u-g(u) 



U 7= 

0, u = 



C\ — sin 



C2M 
M 2 + l 



where < c x = sin f < c 2 < ^4 < 0.02. 



2 ^ ^ ^ 4+3V2 

First of all, notice that for any £ in a small interval centered in zero, the 
2 , n , ^ . Ii^l 

z = 16 n i ^2 
9 



algebraic curve E £ = 1 + ip 2 — (1 + e 
and with the smooth peek %p e + = ^ (1 + e 



is concave in quadrant I 
Notice further that the energy 
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associated to the flow of / can be estimated by 

E(^fi) = y + y /(«)** 

32 



and respectively 

E(ip,S)>'—-^-(l + Cl ) 



So, lying in between the curves E Cl = and E Cl _ C2 = 0, the algebraic curve 
E = will obey the "+/— " splitting condition. 

Second, let us find the zeros of /. Obviously, they are fixed points of g. 
Besides the null solution, the algebraic equation g{u) = u reads as 

2 

1 + c\ — sin — = vH, u 4^ 0. 

ir + 1 

In (0, +oo), the left-hand member of the equation is a decreasing function 
while the right-hand member is increasing (strictly). In consequence, the 
equation can have at most one solution here. Since C\ = sin y, the solution 
reads as Uq — 1. Similarly, we find the third fixed point u± = —Uq. 

Third, to find the quantity 77(a) that describes the growth of the vorticity 
function / around a, notice that 

l/HI < |«| + V\u\(l + Ci + C 2 ) < \ U \ + f 1 + 



2 



and so the inequality |/(ix)| < 77a will reduce to (l + < I ±) m xh e 
function i] i— >■ ^(77) = - ^ +! / is increasing in (3, |]. By noticing the elemen- 
tary inequality § (y^i - l) > > we deduce that 1 + < y/§| ~ 1 < 



a (y), meaning that 77(a) 



10 
3 • 



Fourth, to find the coefficient L(a) of the locally Lipschitzian (outside 
{0}) function g, we have the estimates 

1 i 
If H| < l + ^=-(l + Cl + c 2 ) + 2c 2 



2 



I + C1 + C2 |^jj 1 

1 H ;= h 2c 2 • 5- ■ ; 

|«| (u 2 + l)l v^TT 
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1 + ci + c 2 , s /2 

< l + 2c 2 + < l + 2c 2 + l + d + ca W- 

g Va 



2a/(!-I) 



< 1 + 2c 2 + ( 1 + ^ ) y/2 < 1 



The latter estimate is equivalent to c 2 < |jj^f ■ 
Fifth, given if> > 0, observe that 

l g(u)du > / y/udu ■ (1 + c\ — sin c 2 ) 
Jo Jo 



^3/2 



2 1+ci— sinc2 



where A 



4 1+sin — siriC2 _ 



2A ' 



Sixth, regarding the quantities c and i> from (Tl4j). observe that 

H 3/2 



< 



< 



[(1 - sinc 2 ) + ci] 



# 2 



i? 2 



i? 2 



1 — sin 



c 2 ^ 2 

V^ 2 + 1 



+ Ci 



< J^-(l+*)< 



i? 2 V 2 



and so c = y, ^ = \. 



References 

[1] A. Constantin, R.S. Johnson, Propagation of very long water waves, 
with vorticity, over variable depth, with applications to tsunamis, Fluid 
Dynam. Res. 40 (2008), 175-211 

[2] A. Constantin, R.S. Johnson, Addendum: Propagation of very long 
water waves, with vorticity, over variable depth, with applications to 
tsunamis, Fluid Dynam. Res. 42 (2010), Art. No. 038901 

[3] A. Constantin, A dynamical systems approach towards isolated vorticity 
regions for tsunami background states, Arch. Rational Mech. Anal., doi: 
10.1007/s00205-010-0347-l 



